Determination of the spin-flip time in ferromagnetic SrRu0 3 from time-resolved Kerr 

measurements 
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We report time-resolved Kerr effect measurements of magnetization dynamics in ferromagnetic SrRuC>3. We 
observe that the demagnetization time slows substantially at temperatures within 15K of the Curie temperature, 
which is ~ 150K. We analyze the data with a phenomenological model that relates the demagnetization time to 
the spin flip time. In agreement with our observations the model yields a demagnetization time that is inversely 
proportional to T-T c . We also make a direct comparison of the spin flip rate and the Gilbert damping coefficient 
showing that their ratio very close to k^Tc, indicating a common origin for these phenomena. 



I: Introduction 

There is increasing interest in controlling magnetism in fer- 
romagnets. Of particular interest are the related questions of 
how quickly and by what mechanism the magnetization can be 
changed by external perturbations. In addition to advancing 
our basic understanding of magnetism, exploring the speed 
with which the magnetic state can be changed is crucial to ap- 
plications such as ultrafast laser-writing techniques. Despite 
its relevance, the time scale and mechanisms underlying de- 
magnetization are not well understood at a microscopic level. 

Before Beaurepaire et al.'s pioneering work on laser-excited 
Ni in 1996, it was thought that spins would take nanoseconds 
to rotate, with demagnetization resulting from the weak inter- 
action of spins with the lattice. The experiments on Ni showed 
that this was not the case and that demagnetization could oc- 
cur on time scales significantly less than 1 ps 1 . Since then 
demagnetization is usually attributed to Elliott- Yafet mecha- 
nism, in which the rate of electron spin flips is proportional 
to the momentum scattering rate. Recently Koopmans et al. 
have demonstrated that electron-phonon or electron-impurity 
scattering can be responsible for the wide range of demag- 
netization time scales observed in different materials 2 . Also 
recently it has been proposed that electron-electron scattering 
should be included as well as a source of Elliott- Yafet spin 
flipping, and consequently, demagnetization 3 . Although Ref. 3 
specifically refers to interband scattering at high energies, it 
is plausible that intraband electron scattering can lead to spin 
memory loss as well. 

Time-resolved magneto-optical Kerr effect (TRMOKE) 
measurements have been demonstrated to be a useful probe of 
ultrafast laser-induced demagnetization 1 . In this paper we re- 
port TRMOKE measurements on thin films of SRO/STO( 111) 
between 5 and 165K. Below about 80 K we observe damped 
ferromagnetic resonance (FMR), from which we determine a 
Gilbert damping parameter consistent with earlier measure- 
ments on SrTi03 with (001) orientation 6 . As the the Curie 
temperature (~ 150-fQ is approached the demagnetization 
time slows significantly, as has been observed in other mag- 
netic systems 4 . The slowing dynamics have been attributed 
to critical slowing down, due to the similarities between the 
temperature dependencies of the demagnetization time and 



the relaxation time 5 . In this paper we develop an analytical 
expression relating the demagnetization time to the spin-flip 
time near the Curie temperature. This provides a new method 
of measuring the spin-flip time, which is essential to under- 
standing the dynamics of laser-induced demagnetization. 
II: Sample Growth and Characterization 
SRO thin films were grown via pulsed laser deposition at 
700°C in 0.3 mbar of oxygen and argon (1:1) on Ti02 termi- 
nated STO(lll)2. A pressed pellet of SRO was used for the 
target material and the energy on the target was kept constant 
at 2.1 J/cm 2 . High-pressure reflection high-energy electron 
diffraction (RHEED) was used to monitor the growth speed 
and crystallinity of the SRO film in situ. RHEED patterns 
and atomic force microscopy imaging confirmed the presence 
of smooth surfaces consisting of atomically flat terraces sep- 
arated by a single unit cell step (2.2 Ain the [111] direction). 
X-ray diffraction indicated fully epitaxial films and x-ray re- 
flectometry was used to verify film thickness. Bulk magneti- 
zation measurements using a SQUID magnetometer indicated 
a Curie temperature, T c , of ^155K. Electrical transport mea- 
surements were performed in the Van der Pauw configuration 
and show the residual resistance ratio to be about 10 for these 
films. 

Ill: Experimental Methods 

In the TRMOKE technique a magnetic sample is excited by 
the absorption of a pump beam, resulting in a change of polar- 
ization angle, A6jf(t), of a time delayed probe beam. The ul- 
trashort pulses from a TkSapph laser are used to achieve sub- 
picosecond time resolution. Near normal incidence, as in this 
experiment, A9a' is proportional to the z component of the 
perturbed magnetization, AM Z . AQk is measured via a bal- 
anced detection scheme. For additional sensitivity, the deriva- 
tive of A8k t) with respect to time is measured by locking into 
the frequency of a small amplitude (^500 fs) fast scanning de- 
lay line in the probe beam path as time is stepped through on 
another delay line. 
rV.l: Experimental Results: Low Temperature 
Fig. 1 shows the time derivative of AQk for an 18.5nm 
SRO/STO(l 1 1) sample for the 16ps following excitation by a 
pump beam, for temperatures between 5 and 85K. Clear fer- 
romagnetic resonance (FMR) oscillations are present, gener- 
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FIG. 1 . Derivative of the change in Kerr rotation as a function of time 
delay following pulsed photoexcitation, for 5<T < 85 K 



5.0 - 



4.0 



in 
t3 

CO 



£ 3.0 

CD 

B 2.0 

I 1-0 
0.0 



- Short time Feature 
■ FMR Oscillations 



20 40 60 80 
Temperature [K] 

FIG. 2. Comparing amplitudes of the short time feature and the fer- 
romagnetic resonance oscillations 



ated by a sudden shift in easy axis direction upon thermal ex- 
citation by a pump beam 6 . This motion is described by the 
Landau-Lifshitz-Gilbert equation with the frequency of oscil- 
lation proportional to the strength of the magnetocrystalline 
anisotropy field, and the damping described by dimension- 
less phenomenological parameter, a. The motion appears as 
a decaying oscillation to TRMOKE. The orientation of the 
anisotropy field, closer to in-plane with the sample surface in 
SRO/STO(lll) than in SRO/STO(001), makes these oscilla- 
tions more prominent when observed with the polar Kerr ge- 
ometry compared to previous measurements. 

Attempting to model the time derivative of A9jf with a 
damped cosine reveals that it cannot be fit by such a function 
for t < 2ps. The feature at short times in Fig. 1 contains 
higher frequency components, whereas the oscillations which 
become clear after 2 ps are at a single frequency. A com- 
parison of the amplitude of the first peak (at t ~.5 ps) with 
the amplitude of the subsequent oscillations (defined as the 
difference between dA8^/dt at the peak at ~3.5 ps and the 
dip at ~5.5 ps), is shown as a function of temperature in Fig. 
2. The constant offset between the two amplitudes indicates 
that dAO#7dt is comprised of a superposition of a tempera- 
ture independent, short-lived component with the longer lived 
damped oscillations. 

Fitting the oscillatory portion of the signal to a damped 
cosine, the temperature dependencies of the amplitude, fre- 
quency, and damping parameter are found, as shown in Fig. 3. 
Comparing these parameters for SRO/STO(l 1 1) to previously 
published work on SRO/STO(001), the frequency is found to 
be somewhat smaller and to change more with temperature. 
Of particular interest is a, which is also smaller in this ori- 
entation of SRO, consistent with the more pronounced FMR 
oscillations. Strikingly, in both orientations there is a dip in 
a around 45K, which is relatively stronger in SRO/STO(l 1 1). 
This further strengthens the link between a and the anoma- 
lous hall conductivity, speculated in that paper, through near 
degeneracies in the band structure 6 . 

IV.2: Experimental Results: High Temeperature 



By taking the time derivative of A8^, the FMR oscilla- 
tions can be followed until they disappear at elevated temper- 
atures, at which point it becomes simpler to look at A<9k than 
its time derivative. Fig. 4 shows AO^ as a function of time 
for the first 38 ps after excitation by the pump laser, for tem- 
peratures between 120K and 165K. A property of a second 
order phase transitions is that the derivative of the order pa- 
rameter diverges near the transition temperature. The peak in 
magnitude of AQk in figure 4, shown in figure 5, can be un- 
derstood as the result of the derivative of magnetization with 
respect to temperature becoming steeper near the Curie tem- 
perature. A strong temperature dependence of the demagneti- 
zation time, tm, is seen, with tm significantly enhanced near 
150K, consistent with previous reports on SRO 4 ' 6 . 

A6jf(t) in Fig. 4, normalized by the largest value of 
A9a'(0 in the first 38 ps, can be fit with the following func- 
tion: 



for t < 
for t > 



Ae mox (t) 
Ae K (t) 

AQ max (t) 
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C - Ae-V™ 



(1) 



where the decay time is tm- The resulting tm is plotted 
as a function of temperature in Fig. 6. Notably, tm increases 
by a factor of 10 from 135K to 150K. Taking the fit value of 
T c = 148. 8K, as will be discussed later, tm is plotted log- 
log as function of reduced temperature, tn — (T c — T)/T c . 
The result looks approximately linear, indicating a power law 
dependence of tm on the reduced temperature. 

V: Discussion of Results: 

Efforts to explain demagnetization have been largely phe- 
nomenological thus far, understandably, given the daunting 
challenge of a full microscopic model. Beaurepaire et al. in- 
troduced the three temperature model (3TM) to describe de- 
magnetization resulting from the interactions of the electron, 
phonon, and spin baths 1 . In 3TM the dynamics are determined 
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FIG. 3. Temperature dependence of (a) Amplitude of oscillations, 
(b) FMR frequency, and, (c) damping parameter 



„ 3.0 

LO 




- 120K 
-125K 

130K 
-135K 

- 140K 
-142.5K 

145K 
-147.5K 



-147.5K 
-150K 

152.5K 
-155K 

157.5K 
-160K 

165K 



10 20 30 

Time [ps] 



40 



FIG. 4. Change in Kerr rotation as a function of time delay following 
pulsed photoexcitation, for 120<T < 165 K 
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FIG. 5. Magnitude of change in Kerr rotation at 38ps as a function 
of temperature 
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FIG. 6. Demagnetization time at high temperature 



by the specific heats of each bath as well as the coupling 
constants between them. Demagnetization can generally be 
described with the appropriate choice of coupling constants, 
providing a guide into the microscopic mechanism. Koop- 
mans et al. also offer a phenomenological description of de- 
magnetization considering three baths, but one that follows 
spin in addition to heat 8 . Spin is treated as a two state system 
with energy levels separated by an exchange gap and Fermi's 
golden rule is used to relate demagnetization to electron scat- 
tering which flips a spin. Equations for coupling constants are 
derived based on parameters such as the density of states of 
electrons, phonons, and spins, the electron-phonon scattering 
rate, and the probability of spin flip at a scattering event. 

In the following we attempt to understand the behavior of 
the demagnetization time near T c with an approach based on 
the two spin state model. A general relationship between the 
laser-induced tm and the spin flip time, t s j, can be derived 
near the transition temperature based on the concept of de- 
tailed balance 9 . In equilibrium, the ratio of the probability of 
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a spin flipping from majority to minority to the reverse of this 
process is the Boltzmann factor, e - A ^/ fe T^ where A ex is the 
exchange energy gap. The time derivative of the number of 
majority and minority electrons can then be written: 



N, 



ma j 



-Nrn.ir), — 



_ N ma j _A CT /it B r 



(2) 



When the sample is thermally excited by a pump beam, the 
electron temperature is increased by 5T e . The rate of change 
of spins is then altered in the following way: 



N„ 



-N„ 



N, 



maj c -A^/k B (T+ST e ) 



(3) 



The demagnetization time is related to the total change in 
spin, AS, from initial to final temperature, where, setting h=l, 
S is defined by: 



S = l/2(N maj - N rnin )/N total (4) 
Assuming that AS, as a function of time, can be written: 



As(t) = [s(T f ) - s(Ti)](i - e -* /TM ; 

the demagnetization time can be written as: 



tm 



AS 

W) 



(5) 



(6) 



where 5(0) is the initial change in the time derivative of the 
spin. 

The total change in spin can be calculated by taking the 
derivative of 5 with respect to T, and multiplying by AT eq , 
the increase in temperature once electrons, phonons, and spins 



have come into thermal equilibrium with each other. 5(T) and 
AS can be written: 



5(T) = - - tank ( -A- 
v ' 2 \2kT 



(7) 



and: 
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where we have relied on the fact that near the transition tem- 
perature, A ex -c kgT and made the approximation that 

ST e <C T for low laser power. In the last equation To 3> 1 
near T c , so only the first term will be considered. 

The quantity 5(0), where 5 = \/2{N my - N m i n )/N tota i, 
can be found by taking the derivative of 5(0) with respect to 
T e , since immediately after excitation the electron tempera- 
ture has increased, but the spin temperature, T, has not. 
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(9) 



N m in ~ h N total- 



Near the Curie temperature N ma j . ,„,,,. - .> 
Using this approximations and equation (O, we find: 



tm 



T r T, 



C T s f 



(10) 



where A' ex is the derivative of A ex with respect to temperature 
and A ex ~ (T c — T) 13 , where j3 is the critical exponent of 
the order parameter. Taking the derivative, we find A' ex ~ 
—(3(T C - T) f3 ~ 1 , and thus can write 



tm 



pT sf 



T r -T 



(ID 



Therefore tm is predicted to scale as 1/ (T c — T) near the 
transition temperature. A fit of T c ~ 148. 8K is found for the 
data in Fig 6. 

Note that detailed balance suggests that the demagnetiza- 
tion time scales as 1/tn near the transition temperature re- 
gardless of the underlying mechanism of the demagnetization. 
Additionally, the critical exponent found is independent of j3. 
It should also be noted that the current situation, where the 
sample has been excited by a laser, is distinct from critical 
behavior as typically considered. In general, divergent time 
scales are linked to divergent length scales, but here excita- 
tions of various length scales are not being excited. Instead 
the length scale is always effectively infinite, having been de- 
termined by the laser spot size. r s f is plotted as a function of 
temperature for the mean field value of f3 = 1/2, which has 
been shown to be suitable for SRO 10 , in Fig. 8. t s j is revealed 
to be approximately 200 fs and nearly constant as a function 
of temperature. 
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FIG. 8. Spin flip time at high temperature 



Previous reports of conductivity in SRO give a scattering 
time of ~20 fs near the transition temperature 11 . A compar- 
ison of the spin flip time with the scattering time implies a 
probability of 0.1 that a scattering events results in a spin flip. 
Though electron-phonon interactions are the most commonly 
considered source of demagnetization, as mentioned previ- 
ously, Eliot Yafet-like electron-electron coulomb scattering 
can also result in demagnetization 3 . This is especially true for 
materials with strong spin orbit coupling, such as SRO. Ad- 
ditionally in SRO the interaction with the crystal field means 
that total spin is not conserved [Goodenough], so every elec- 
tron interaction can perturb the spin state. 

Having found a relationship between the demagnetization 
time and the spin flip time we would like to explore the rela- 
tionship between these parameters and the damping param- 
eter, a. Intuitively, the damping parameter should be pro- 
portional to the spin flip scattering rate, or inversely propor- 
tional to the spin flip scattering time: a ~ 1/t s /. Elliot- Yafet 
type scattering dissipates energy from motion described by 
the LLG equation by disrupting the coherent, collective pre- 
cession of spins. Spins that have had their angular momen- 
tum changed through electron collisions must be pulled back 
into the precession through the exchange interaction, repre- 



senting a transfer of energy away from the precessional mo- 
tion. These collision-mediated spin-orbit coupling effects are 
thought to be the primary source of Gilbert-type damping in 
ferromagnets 12 . Again, this should be particularly true in a 
ferromagnet with strong spin orbit coupling. 

Combining the spin flip time and the damping parameter 
with Planck's constant reveals an energy scale, £, given by 
the condition that: 



1 

a 



Tsf 



(12) 



Noting that the values for a and t s j found in figures 3 and 
7, respectively, are approximately constant as a function of 
temperature, this energy scale for SRO is ~7 meV. The fun- 
damental energy scales applicable to the magnetic system in 
SRO are the Fermi energy, the exchange energy, and the criti- 
cal temperature, the last two of which are interdependent. The 
Fermi energy is orders of magnitude larger than 7 meV, but 
the energy associated with the critical temperature, ksT c ~ 
13 meV, is of the same order. This suggests an underlying 
connection between the critical temperature (and thus the ex- 
change energy), Gilbert damping, and spin flip scattering. 

A relationship similar to equation (fT2l) has been found pre- 
viously between tm (rather than r s f) and a by Koopmans et 
al. at low temperature: 



1 h 1 



tm = t 



4 ksT c a 



(13) 



Applying this equation to SRO at 5K yields r m ~ 30fs, 
which is unphysical since it is below the total scattering rate 
of ~100fs at low temperature 11 . Whether the fundamental re- 
lationship is between transition temperature and the demagne- 
tization time or the spin-flip scattering time remains a question 
for a microscopic model to resolve. 
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